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SIGNAL PROCESSING APPARATUS AND METHOD 
Field of the Invention 

The present invention relates to signal processing apparatus and methods for 
representing the effects of radio communications channels. 
5 Background of the Invention 

Data is communicated using radio signals by modulating a radio frequency 
signal with the data in some way. The radio frequency signal propagates from a 
transmitting antenna and is detected by a receiving antenna. An estimation of the 
communicated data is then recovered by de-modulating the received radio signal. 

10 In order to develop techniques for communicating data using radio frequency 

signals, models which represent the effects on the data during communication from the 
transmitter to the receiver antenna have been developed. For example, for mobile 
radio telephone systems, propagation and channel modelling studies have been 
performed for radio frequency bands utilised by the Global System for Mobiles (GSM) 

15 and the Universal Mobile Telecommunications System (UMTS) for 3 rd Generation 
mobile (3G). To develop such models channel sounding measurements were made and 
in combination with theoretical considerations, statistical channel models were 
developed. 

Radio communications channels associated with mobile radio systems such as 
20 GSM and UMTS are characterised by including a number of multi-path rays with 
distinct propagation delays. Each path represents the propagation of a version of a 
radio signal from a transmitter to a receiver. It is known that the amplitude of each 
version of the radio signal received via each of the paths can be represented as zero 
mean complex gaussian random variable. At the receiver the multi-path versions 
25 combine to form a composite signal. If the difference between the respective delays 
for each path exceeds a symbol period of the communicated data, then inter-symbol 
interference occurs, which must be corrected in order for data to be communicated. 
Hence, receiver techniques are arranged to recover data in the presence of multi-path 
propagation. In order to simulate a multi-path radio communications channel, a 
30 complex gaussian random variable is generated for each path and for each sample of 
the simulated radio signal. Therefore, the simulation of such multi-path channels can 



P016130EP 2 

represent a substantial computational task, requiring a relatively long time to generate 
simulation results representing the performance of a radio communications receiver. 
In particular, but not exclusively, a time taken to simulate spread spectrum 
communications can be prohibitively long. This is because a spread spectrum signal is 
5 formed by arranging for the data to be communicated to modulate a spreading code. 
The spreading code includes a much greater number of bits, referred to as chips, than 
the data being communicated. Thus for each chip a multi-path fading channel must be 
simulated, representing a greater computational task than that for simulating a base 
band data symbol, where no spreading is modelled. For this reason it is desirable to 

10 reduce the number of computations by simplifying the simulation of a radio 
communications channel. 

Examples of communications techniques, which might be simulated, include 
channel coding techniques such as forward error correction encoding and decoding. In 
order to evaluate such techniques for a radio communications channel it is desirable to 

15 represent the time correlation effects of the radio communications channels. This is 
because radio communications channels are characterised by fading effects, where 
multiple rays of the radio signal combine destructively during some instances and 
constructively at other instants, which can have an effect of introducing bursts of 
errors of the communicated data. Therefore the performance of the radio channel is 

20 time correlated. 

A technical article by Stefan A. Fechtel entitled "A Novel Approach to 
Modeling and Efficient Simulation of Frequency-Selective Fading Radio Channels," 
published in the IEEE Journal on Selected Areas in Communications, Vol. 11, No. 3, 
April 1993, discloses a two step channel orthogonalisation technique. The 

25 orthogonalisation technique only provides a facility for obtaining a bounded 
performance for channel coding techniques, whereas a more exact performance 
representing the time correlation effects of the radio channel is preferable. 
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Summary of Invention 

According to the present invention there is provided a signal processing 
apparatus operable to represent the effects on a received signal of a radio 
5 communications channel having L paths. Each path has an average attenuation and a 
pre-determined respective delay. The received signal includes a combination of 
correlated components determined from an effect of pulse shaping filters on the 
received signal, each correlated component having a correlation coefficient 
representing a correlation of the received signal component with respect to each of the 

10 other components. The signal processing apparatus comprises a plurality of signal 
simulators, each simulator generating a signal component value proportional to a 
complex zero mean gaussian random variable having a pre-determined variance. The 
signal processing apparatus includes a summer operable to sum the signal component 
values produced from each signal simulator, to form a representation of the signal 

15 received via the radio communications channel. The variance of each of the signal 
simulators is pre-determined by calculating the eigen values of a matrix formed from 
the correlation coefficients and from a channel correlation matrix which includes the 
average attenuation of each of the L paths. 

Embodiments of the present invention can provide simulators for radio 

20 communications channels, which can represent the radio communications channel with 
a substantially reduced complexity with respect to known simulators. The reduction in 
complexity is achieved by performing a transformation of a conventional 
representation of a radio communications channel having L paths and L correlation 
coefficients. As a result of the transformation the radio communications channel may 

25 not only be simulated with substantially reduced complexity but also removes a 
requirement for representing the correlation between correlated components from 
which the received signal is formed. Furthermore in some embodiments, the number 
of signal simulators may be less than the number of paths L of the radio 
communications channel being simulated. This is because after the transformation has 

30 been performed, the effect of the signal simulators with longer delays is relatively 
small compared with signal simulators with shorter delays and so ignoring these terms 
can provide a good approximation to the simulated channel. Therefore in such 
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embodiments, the number of signal simulators is determined from the number of eigen 
values above a pre-determined threshold, each eigen value above the threshold 
forming the pre-determined variance for a corresponding signal simulator. 

In some embodiments, the correlation coefficients, representing a correlation 
5 between each correlated component with respect to each of the other components, may 
represent a correlation between the output of each correlator of a rake receiver. For 
Code Division Muliple Access (CDMA) communications a rake receiver is often used 
to recover data from a spread spectrum signal. Each correlator of the rake receiver 
correlates the received signal with respect to a spreading code at a pre-determined 

1 0 delay. The output of each rake correlator may be correlated with respect to the outputs 
of the other correlators, the correlation being represented by the correlation 
coefficients. However, embodiments of the present invention are not limited to 
simulating CDMA communications. In other embodiments the separate components 
of the communicated signal may be formed Multiple Inputs Multiple Outputs (MIMO) 

1 5 communications or Time Division Multiple Access (TDM A). 

Various further aspects and features of the present invention are defined in the 
appended claims and include a method of representing the effects of a radio 
communications channel, a method of simulating a radio communications channel and 
a computer program. 

20 
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Brief Description of the Drawings 

Embodiments of the present invention will now be described by way of 
example only with reference to the accompanying drawings, where like parts are 
provided with corresponding reference numerals, and in which: 
5 Figure 1 is a schematic block diagram of a transmitter and receiver chain of 

parts which are typically involved in simulating data communications via a radio 
communications channel; 

Figure 2 is a schematic block diagram of a multi-path correlation receiver 
otherwise known as a rake receiver; 
10 Figure 3 is a schematic block diagram of the transmitter and receiver chain of 

Figure 1 adapted to include a simplification of parts representing the radio 
communications channel; 

Figure 4 is a schematic block diagram of a signal processing apparatus which 
simulates the parts of the communications channel illustrated in Figure 3; 
15 Figure 5 provides a graphical representation of simulation results for a 

"Pedestrian A" radio communications channel with a mobile motion speed of 3km/h, 
for the simplified equivalent channel and the conventional fully simulated channel; 

Figure 6 provides a graphical representation of simulation results for a 
"Pedestrian A" radio communications channel with a mobile motion speed of 
20 120km/h, for the simplified equivalent channel and the conventional fully simulated 
channel; 

Figure 7 is a tabular representation of the components making up a 3GPP 
"Typical Urban" channel, before and after transformation; and 

Figure 8 provides a graphical representation of multi-path channel components 
25 before and after transformation according to the simplification; 

Figure 9 is a flow diagram illustrating the operations involved in simulating the 
communications channel illustrated in Figures 3 and 4. 



30 
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Description of Preferred Embodiments 

Schematic block diagrams in Figures 1, 2 and 3 provide an illustration of a 
transformation providing a simplification of a simulated radio communications 
channel according to an embodiment of the present invention. In Figure 1 a data 
5 source 1 generates data in digital form, which is fed to an encoder 2. The encoder 2 
could apply any form of coding to the data generated from the source 1 , such as error 
correction encoding. A spread spectrum encoder 3 then receives the encoded data, and 
generates a spread spectrum signal by modulating a spreading code with the encoded 
data. Typically, the spreading code is formed from a pseudo-random bit sequence 

10 having a non-repetitive period of one thousand or more bits, each bit of the modulated 
spreading code being known as a chip. The spread spectrum signal is then received by 
a modulator 4. The modulator modulates a base band carrier signal with the chips of 
the spread spectrum signal. Typically this might include some form of phase or 
amplitude modulation such a QPSK or QAM in which the phase and/or amplitude of 

1 5 the base band carrier signal is modulated with the encoded data. 

The signal from the modulator 4 is fed to a pulse shaping filter 6 which 
provides a bandwidth limiting effect which would be applied to the radio signal to be 
communicated before transmission. As those acquainted with Nyquist bandwidth 

20 filtering will appreciate, the pulse-shaping filter is typically in the form of a root raised 
cosine response so that by employing a corresponding root raised cosine filter at the 
receiver, an overall raised cosine response is produced. Accordingly by sampling the 
received data at the symbol time (which for a spread spectrum signal is the chip 
period), the effects of adjacent chips will be substantially zero at the sampling instant, 

25 thereby minimising inter-symbol interference which otherwise might be introduced by 
the band limiting filters. 

For an equivalent transmitter to that represented in Figure 1, the base band 
modulated signal would be up converted to a radio frequency signal, amplified and 
transmitted from an appropriate antenna. However, as those acquainted with the 

30 simulation of digital communications systems will appreciate, the effects of the radio 
communications channel can be modelled by applying equivalent effects to the base 
band signal produced by the modulator 4. 
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Many radio communications channels can be represented as a plurality of L 
discrete paths, each path providing an independently fading version of the radio signal 
being represented by a zero mean complex gaussian process also known as a Rayleigh 
fading process.. Therefore as shown in Figure 1, the simulated transmitted signal 
5 passes from the pulse-shaping transmitter filter 6 to an anti-aliasing receiver filter 8, 
via an L-path Rayleigh fading channel 10 and a noise generator 12. Whilst the L-path 
Rayleigh fading channel 10 models the effects of a multi-path radio channel, the 
effects of noise which is introduced during transmission and is present at the receiver 
antenna (due to thermal noise) is provided by the noise generator 12. If u is the signal 
10 before channel 10, then the signal after the noise generator 12 is 

L-\ 

v(t) = Y d h i u(t-r i ) + AWGN 

with hj the complex channel coefficient and t\ the delay for path i. AWGN stands for 
Additive White Gausssian Noise. 

After passing through the receiver filter 8, the received signal is fed to a de- 

15 modulator 14, which produces an estimate of the spread spectrum encoded data 
produced by the spread spectrum encoder 3, by performing a reverse mapping of the 
data from the complex base band signal as was applied to the modulator 4. The 
received spread spectrum encoded data is then fed to a multi-path correlation receiver 
16. As will be explained shortly with reference to Figure 2, the multi-path correlation 

20 receiver 16, recovers an estimate of the base band encoded data, by correlating the 
received demodulated data with a reproduced version of the spreading code which was 
used to form the transmitted signal by the spread spectrum encoder. The de-spread 
base band encoded data is then fed to a decoder 18, which forms a reverse of the 
encoding process to estimate the base band date. The estimation of the base band data 

25 may include correcting errors in the estimation of the received data, if the data has 
been encoded using an error correction code. The estimated base band data is then fed 
to a sink 20. 

Figure 2 provides a more detailed representation of the multi-path correlation 
receiver, otherwise known as a rake receiver. The rake receiver 16 includes a plurality 
30 of L correlators 30, each of the correlators 30 receives on a first input 32 a version of 
the spreading code C n (t) which was used to spread the received data bearing spread 
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spectrum signal in the spread spectrum encoder 3. To a second input 34 of the 
correlators 30, the received sampled signal is fed from the de-modulator. The samples 
of the received signal are multiplied by the corresponding samples of the spreading 
code at the chip rate. The result of being multiplied by the spreading code is summed 
5 by summers 36 to form, for each spread spectrum symbol, a de-spread sample, which 
is summed by a further summer 38. The further summer 38 forms a sample of the base 
band encoded data which is de-spread and has the effect of substantially reducing the 
inter symbol interference introduced by the multi-path fading channel. 

The spreading code is reproduced within the rake receiver 16. Each version of 

10 the spreading code fed on the respective inputs 32 for each correlator 30 is delayed 
with respect to the first correlator, by an amount which represents a likely temporal 
position of one of the paths of the L path channel. Accordingly, the delays of each of 
the spreading codes which are represented as corresponding shifts of the spreading 
code sequence, aim to reflect the distribution of energy introduced by the multi-path 

15 channel. 

From the representation of the spread spectrum communications channel 
according to the example of Figure 1 , a base band data symbols s communicated from 
the output of the encoder 2, to the output r of the multi-path correlation receiver 18, 
may be represented by the following equation: 

20 

r = J JSZ^>» +ISI + AWGN (1) 



Each of the L paths of the simulated multi-path channel have parameters 
(4 > r , )fl < i <h where X t is the average attenuation of path / with a delay The 
25 received symbol r is determined for a transmitted symbol s from equation (1) in which 
X t is a complex zero-mean gaussian random variable of variance X t and py is the 

correlation between the correlators i and j of the rake receiver, the number of 
correlators in the rake receiver being L. For the example of root raised cosine filters, 
the correlation coefficients can be derived from the delays through the auto-correlation 
30 function of the raised cosine filter, as expressed below: 
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sin(;rA (y ) cos(7rfiA ij ) 
^(1-4^4) 



(2) 



where 




and for the example of UMTS, the parameters for these 



equations are: 



/? = 0.22 
V T = 3MMHz 



10 



15 



20 



As mentioned above, in order to simulate the effects of the multi-path fading 
channel, the Z-path multi-path fading channel 10, requires the generation of a complex 
gaussian sample per chip of the spread spectrum signal. As will be appreciated, this 
represents a substantial computational load, in terms of the number of computations 
per second. Accordingly, a time taken to simulate the communication of a base band 
symbol including spreading and de-spreading, with multi-path fading channel at the 
chip rate may be large and for some simulations prohibitive. 
Transformation of L-path Multi-path Fading Channel 

Embodiments of the present invention provide a facility for reducing the 
number of computations required per second required to simulate a multi-path fading 
channel. The reduction in complexity according to an embodiment of the invention 
comprises introducing a transformation of a matrix of values representing the double 
summation expressed as equation (1) above. The transformation thereby provides a 
simplification of the channel model. 

The transformation comprises computing the eigen values (Mi)o<i<i°^ a matrix 
formed according to the double summation of equation (1), which is [( p.. ,) 0 ^ iJ<L \K 
where K is the correlation matrix of the channel (channel correlation matrix). In the 
usual case of independent paths, it becomes [(p u ) 0 ^ J<L J- ^g[(^ ) 0 <,,y</J- According 
to the simplification provided by the transformation, the schematic block diagram of 
Figure 1, becomes that shown in Figure 3, where the output of the rake receiver 16' 
according to the simplified parts shown within a box 50, becomes that expressed by 
the equation below: 
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JAM 2 
+ 1 SI + A WGN (3) 

Where in equation (3), Y( is a complex zero-mean gaussian random variable of 

variance /// and there is no correlation component. The main embodiment of the 

5 present invention is therefore to replace the channel definition set (k\, xj) by (^ij, iT c ). 
For any initial channel model, we find a equivalent one which is adapted to the rate of 
the transmission we consider. It can be CDMA, TDMA (replacing the chips by 
symbols) or MIMO. Furthermore, according to equation (3), and as explained below, 
the number of signal simulators N may be less than or equal to the number of paths L. 

10 The radio communications channel represented by the schematic block diagram 

shown in Figure 3, corresponds to that shown in Figure 1 and so only the differences 
between Figures 1 and 3 will be explained. According to the transformation expressed 
by equation (3), the multi-path Rayleigh fading channel is transformed, to form an 
equaivalent channel 10'. As a result of the transformation all parts within the box 50 

15 may be represented by equation (3). A diagram providing an implementation of a 
signal processing apparatus for representing the parts in box 50 is shown in Figure 4. 

In Figure 4, a plurality of N signal simulators 60 generate signal component 
values Yf in accordance with zero mean complex random gaussian processes having a 

variance /up The component values are fed to data processors 62 which form the 

20 squared magnitude |7| 2 of the signal component values which are fed to a first input 

64 of multipliers 66. To a second input 68 the input data signal s is fed. The data 
signal is scaled by each of the signal component values by the multipliers 66 and fed to 
a summer 70, which forms the representation of the received signal r. Effectively, the 
signal component values generated from each of the complex gaussian process form 
25 weighting factors for weighting the input data signal. 

Figure 5 and 6 provide a comparison of results produced by a conventional full 
simulation of the multi-path fading channel with respect to results produced using the 
simplified equivalent channel. The results are for a "Pedestrian A" channel simulated 
with a data rate of 12.2 kbps and a motion speed of 3 km/h for Figure 5 and a motion 
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speed of 120 km/h for Figure 6. Bit Error Rate (BER usual) and Block Error Rate 
(BLER usual) results for the usual simulation channel are illustrated with respect to Bit 
Error Rate (BER simplified) and Block Error Rate results (BER simplified) for the 
simplified equivalent channel. As can be observed there is a close correspondence 
5 between the results. 

Truncation of the Number of Simulated Paths 

Equation (3) formed using the transformation according to an embodiment of 
the invention is less complex to compute than equation (1), amongst other reasons 
because there is no double summation. Furthermore, typically as a result of the 
10 transformation, a further simplification can be introduced. This is because the later 
values for ft (those closer to L) are typically low and have a small influence on the 

resulting calculation of the effects of the channel on the received symbol r, than the 
earlier components closer to zero delay. Therefore, even if it is not an exact result, a 
good approximation results when considering only the ft values over a certain pre- 

15 determined threshold. As a result, a further simplification of the simulation of the 
multi-path channel may be produced, by not including signal simulators generating 
complex zero mean gaussian having eigen values which are below a pre-determined 
threshold. Therefore in such embodiments, the number of signal simulators N is 
determined from the number of eigen values above the pre-determined threshold, each 

20 eigen value above the threshold forming the pre-determined variance for a 
corresponding signal simulator. 

An example of this further simplification is provided in table 1 shown in Figure 
7. In the first column labelled "before", the average attenuation with respect to delay 
for a "Typical Urban" channel adopted by the 3 GPP is shown. In the second column 

25 labelled "after", the variance of the signal simulators used in the simplified equivalent 
channel is shown after transformation at the corresponding delays. Figure 8 provides a 
graphical representation of the values shown in the table of Figure 7. As can be seen 
from Figures 7 and 8, the number of signal components simulated can be halved using 
the simplified equivalent channel model, by ignoring the signal components which 

30 have a variance ft which is less than 20dB and therefore do not contribute significantly 
to the simulation results. 
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Simple Worked Example 



10 



For a multi-path channel characterised by three paths, with the following 
parameters: 

At r 0 =(U 0 =0.943; 

Atr, =110/15,^ =0.309; 
Atr 2 = 190ns, X 2 =0.103; 

For correlation coefficients describing the reception of a radio signal via the 
multi-path channel, as follows: 

1 0.725 0.318" 
K/>,)o^ 2 ]= 0.725 1 0.848 
0.318 0.848 1 



Then 



0.943 0.224 0.0327' 
0.684 0.309 0.0873 
0.3 0.262 0.103 



for which ju t = 



1.153 
0.197 
0.006 



15 Proof of Equivalents of the Simplified Simulated Channel 

Perfect channel estimation is assumed and inter-symbol interference is 
neglected. Some notations are specific to this proof. X represents the first path and Y 
the second one instead of X 0 and Xi in the rest of the document. The subscripts are 
used for the time scale (symbol k is written Sk ). After the transformation, X' and Y' 
20 are used instead of Y 0 and Y\ in the rest of the document. In the case of two paths, the 
usual channel definition gives (X h 0) for the first path and {X 2 , r) for the second one. 
The correlation coefficient between the two paths is derived from x and the 
characteristics of the pulse shaping filter. 

The symbol k at the output of the Rake receiver can be written 
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r k =s k Je~ k +AWGN 



The energy of this symbol is 



e k =\X k \ 2 +\Y k \ 2 +2p i2 X k Y k 



where Xk and Y k are two independent zero-mean Gaussian processes with 
variance X\ and \i. Diagonalizing the product of the quadratic form and the two paths 
covariance matrix, leads to the eigenvalues \i\ and We consider the energy 



where X\ and Y'k are two independent zero-mean Gaussian processes with 
variance jlai and The associated symbol is 



We want to prove that the two processes e and e' are identical which implies in 
particular r k = r\ and the equality of the error probabilities after decoding. In order to 
achieve this, it is enough to prove that the joint densities p(s) and p(e') with 8 = (si ... 
s n ) and e' = (e'i ... e' n ) are equal for all n. The following demonstration assumes n = 2 
without loss of generality. 

The characteristic function of the symbol energy can be written 



s\=\X\\ 2 + \Y\\ 2 



r' k =s k fi\+AWGN 



E[exp(jv l e l +jy 2 £ 2 )] 
exp[jv* e - -(X * K~ l X + Y* K;'Y)] 



J 4x 2 JdetK x det K Y 

\^[u*{G v -K^)u]du 



dXdY 



with 
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10 



X = {X x ,X 2 ),y x =E[Xt) = E[X 2 2 ],Y = (Y i ,Y 2 ),y r = E[Y?} = E[Y 2 2 }, 



Yx 

\^E[X,X 2 ] 
u = (X l ,Y l ,X 2 ,Y 2 ), 



E[X,X 2 ^ 



Yx J 



Yy 
IW 2 1 



E[YJ 2 i\ 

Yy 




_ E[X X X 2 ] _ E[Y,Y 2 ] 

Yx Yy 



( 



\ 



, 0 jv 2 Gj 
The same method leads to 



a: 



v e iy) = 



J X K 
K 



1 



det(/-A^) 



with 



Kxr = 



Jv 2 J { D jv 2 D j 

where D is the diagonalization of GK so that GK = P*DP. 
Furthermore, a simple calculation shows that 



G v K xy = 



( p* 


°] 




(P 


°1 


v0 


p * 

r J 


KxY 


,0 


Pj 



therefore 



det(/-G |/ ^) = det(/-A KVy ) 



15 



which implies 

Finally the two processes are the same. The simplified one obtained with the 
invention is equivalent to the usual one. So the error probabilities after decoding are 
the same. It is easy to generalize this proof to the case of more paths. In the case of 
imperfect channel estimation, a variant of this proof leads to the same result. 
Other Applications 



One advantage provided by embodiments of the present invention is that the 
simplified equivalent channel model can be used to provide a symbol level Markov 
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chain modelling the parts within the box 50 of Figure 3, including spreading, 
transmission via the communications channel and de-spreading. A Markov channel 
model can represent a radio communications channel using a plurality of channel 
states representing the states of the radio communications channel. A transition 
5 between states is determined according to transition probabilities. It is the transition 
probabilities, which are determined from the effects of the radio communications 
channel represented, for example, by the signal processing apparatus illustrated in 
Figure 4. An explanation of a process for generating the transition probabilities from 
the results of a simulated channel is provided in more detail in a technical article by 
10 A. Saadani, and P. Tortelier, entitled "A First Order Markov Chain Based Model for 
Flat Fading Channel," published at the International Symposium on Personal, Indoor 
and Mobile Radio Communications, PIMRC'2002. 
Summary of Operation 

In summary a flow diagram presenting in Figure 9 summarises a process of 

15 modelling the effects of a radio communications channel. Starting with step SI, 
channel sounding is performed to identify a number of paths L via which a signal may 
be received from the radio communications channel. The channel sounding process 
also determines at step S2 an average attenuation and a pre-determined delay with 
respect to a first of the paths of a communicated radio signal for each of the paths. 

20 Generally, the steps SI and S2, which characterise the radio communications channel, 
are performed separately to the channel modelling and simplification of the channel 
model according to embodiments of the present invention. Accordingly, in some 
embodiments these steps are omitted or have already been performed for channels for 
which the characterising parameters of the number of paths L, the delay and the 

25 average attenuation have already been established, and so may be omitted. This is 
represented in Figure 9 by a dashed line between steps S2 and S4. 

The steps involved in representing the effects of the radio communications 
channel are therefore represented in Figure 8 as follows: 

S4: A plurality of correlation coefficients representing a correlation between 

30 each of a plurality of components of the received signal are determined in accordance 
with a communications technique which is used to communicate the received signal. 
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For example, the communications technique may be CDMA, in which the coefficients 
represent the correlated output of each finger of a rake receiver. 

S6: A matrix is formed from the correlation coefficients and the average 
attenuation of each of the paths ([(^) 0 ^ ; <J' £'tfg[Wo*.y< J)- 
5 S8: For each of the paths of the radio channel, a variance of a complex zero 

mean gaussian random process is calculated from the eigen values of the formed 
matrix. 

S10: For each simulated signal component N (where N is less than or equal to 
the number of paths Z,), a value from a complex zero mean gaussian random variable 
10 having the variance calculated in S8. 

SI 2: Form the squared magnitude of each of the signal component values 
calculated in S10. 

SI 4: The signal component values are summed to produce, for each path, a 
representation of a signal received via the radio communications channel. 

1 5 As will be appreciated, once the variances of the zero mean complex gaussian 

random process have been determined in S8 from the eigen values of the matrix 
formed in S6, then these steps are no longer necessary when the channel is being 
simulated. Hence in Figure 9 steps S8 and S10 are connected by a dashed line. 
Accordingly, in some embodiments only steps S10, S12 and S14 are performed in 

20 order to represent the effects of the radio communications channel being modelled. 

Various modifications may be made to the embodiments described above 
without departing from the scope of the present invention. For example, the present 
invention finds application in the simulation of receiver diversity and equalisation 
techniques and can be applied to various communications techniques such as TDMA 

25 and MIMO. Various further aspects and features of the present invention are defined 
in the appended claims. 



